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We  consider  first  a  dielectric  medium  of  identical  two-state  atoms  coupled  by  the  radiation  field 
to  an  initially  excited  atom  outside  the  dielectric.  From  the  Schrodinger  equation  follows  a  delay- 
differential  equation  describing  how  the  atom  interacts  with  the  dielectric  by  virtual  photon  ex¬ 
changes.  In  the  macroscopic  limit  of  a  continuous  distribution  of  atoms  in  the  dielectric,  we  derive  a 
simpler  delay-differential  equation  in  which  a  Fresnel  reflection  coefficient  appears.  We  apply  our 
results  to  a  model  of  an  atom  in  a  multimode  Fabry-Perot  resonator,  and  obtaip  a  general  delay- 
differential  equation^  tor  the  probability  amplitude  of  the  initially  excited  state.  This  equation 
predicts  well-known  Rabi  oscillations  when  the  round-trip  photon  propagation  time  is  negligible 
compared  with  the  inverse  of  the  Rabi  frequency  and  the  mirrors  are  highly  reflective.  For  low  mir¬ 
ror  reflectivities  we  recover  Purcell's  prediction  that  the  emission  rate  is  enhanced  by  the  cavity  Q 
factor.  When  the  photon  bounce  time  is  large  compared  with  the  inverse  Rabi  frequency,  Rabi  os¬ 
cillations  do  not  occur.  We  discuss  the  Ewald-Oseen  extinction  theorem  from  the  standpoint  of 
quantum  mechanics. 
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I.  INTRODUCTION 

In  recent  years  there  has  been  considerable  experimental 
and  theoretical  interest  in  the  effects  of  cavity  walls  on 
atomic  absorption  and  emission  processes.  Such  effects 
include  the  inhibited  absorption  of  blackbody  radiation1 
and  the  enhancement2  and  suppression’  of  spontaneous 
emission  rates.  The  theory  of  such  cavity  effects  seems, 
by  and  large,  well  understood.4 

Of  related  theoretical  interest  is  the  Jaynes-Cummings 
model  in  which  a  two-state  atom  interacts  with  a  single 
mode  of  the  electromagnetic  field  in  the  dipole  and 
rotating-wave  approximations.  Jaynes  and  Cummings5 
showed  that  the  same  sort  of  “Rabi  oscillations”  are 
predicted  regardless  of  whether  the  field  is  treated  classi¬ 
cally  or  quantum  mechanically.  Their  paper  was  perhaps 
the  first  detailed  exposition  of  what  is  now  called  the 
“dressed-atom"  approach  to  resonant  atom-field  interac¬ 
tions.  More  recently  Eberly  et  al.b  have  stimulated 
renewed  interest  in  the  Jaynes-Cummings  model  by 
describing  certain  “collapse  and  revival”  effects  that  have 
now  apparently  been  observed  by  Walther  et  a l.1 

Classical  analogs  of  such  effects  follow  from  considera¬ 
tion  of  an  oscillating  dipole  inside  a  cavity.  The  cavity 
walls  reflect  radiation  back  to  the  dipole,  altering  its  oscil¬ 
lations  and  therefore  its  rate  of  radiation.  One  question 
we  address  in  this  paper  is  how  such  multiple  reflections 
of  radiation  off  the  cavity  walls  are  described  quantum 
mechanically.  For  instance,  if  at  time  t  =0  an  atom  is 
suddenly  injected  into  a  cavity,  it  does  not  “know”  it  is  in 
a  cavity  until  a  time  2 d/c,  where  d  is  the  distance  to  a 
cavity  wall.  Such  an  effect  cannot  be  described  within  the 
Jaynes-Cummings  model  in  which  one  begins  a  priori 
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with  a  single-mode  model,  for  retardation  involves  many 
field  modes  in  an  essential  way.  In  our  approach,  there¬ 
fore,  the  single-mode  (Jaynes-Cummings)  results  are  de¬ 
rived  without  the  a  priori  assumption  of  a  single-mode  in¬ 
teraction. 

Another  question  concerns  cavity  damping.  Purcell8  in 
1946  argued  that  for  a  lossy  cavity  the  spontaneous  emis¬ 
sion  rate  should  be  increased  by  a  factor  Q ,  the  cavity 
quality  factor.  Sachdev9  has  considered  this  problem  in 
the  single-mode  context,  and  has  shown  that  Purcell's  pre¬ 
diction  is  justified  in  the  case  of  an  overdamped  cavity;  in 
the  underdamped  case  the  Rabi  oscillations  are  recovered, 
but  they  are  damped  by  the  factor  e~r',  where  y  is  the 
field  damping  rate.10  A  feature  of  our  approach  here  is 
that  the  field  loss  is  not  distributed  but  is  lumped  at  the 
mirrors  as  a  consequence  of  imperfect  reflectivity. 

Since  we  intend  our  approach  to  be  fully  quantum 
mechanical,  we  wish  to  show  how  the  reflection  coeffi¬ 
cient  follows  from  the  Schrodinger  equation  describing 
the  coupling  of  the  atoms  of  the  (dielectric)  mirror  to  the 
field.  This  we  do  in  Sec.  II.  In  Sec.  Ill  we  apply  the  re¬ 
sults  to  an  atom  in  a  Fabry-Perot  resonator  of  length  L 
and  mirror  reflectivities  R.  We  derive  a  delay-differential 
equation  describing  an  atom  in  a  multimode,  lossy  cavity. 
When  flr«l,  where  O  is  the  Rabi  frequency  and 
T  =  2 L/c  is  the  photon  bounce  time,  we  recover  known 
single-mode  results.  When  fl7'»l,  on  the  other  hand, 
the  initially  excited  atomic  state  decays  exponentially  with 
no  Rabi  oscillations.  In  Sec.  IV  we  discuss  the  Ewald- 
Oseen  extinction  theorem  and  summarize  our  results.  Our 
goal  in  this  paper  is  mainly  to  understand  how  dielectric 
mirrors  may  be  described  in  a  fully  quantum-mechanical 
way. 
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II.  ATOM  NEAR  A  PLANE  DIELECTRIC 
INTERFACE:  MICROSCOPIC  THEORY 


We  consider  first  the  situation  illustrated  in  Fig.  1.  At 
xt)  is  an  initially  excited  two-state  atom  with  transition 
frequency  w().  It  is  near  a  dielectric  slab  of  N  ground- 
state  atoms  per  unit  volume,  each  of  which  has  transition 
frequency  «  =  For  simplicity  we  assume  that  all  the 
atoms  have  the  same  (real)  transition  dipole  moment  fi. 

Let  bit)  be  the  probability  amplitude  for  the  state  in 
which  the  atom  at  x0  is  excited,  all  other  atoms  are  in 
their  ground  states,  and  the  field  is  in  its  vacuum  state  of 
no  photons.  Similarly  let  bj(t)  be  the  amplitude  for  the 
state  in  which  the  atom  at  x;  (^=x0)  is  excited,  all  other 
atoms  are  in  their  ground  states,  and  no  photons  in  the 
field.  Finally  let  an(t )  be  the  amplitude  for  the  state  in 
which  all  the  atoms  are  in  their  ground  states  and  the 
field  contains  one  photon  in  the  nth  mode.  Then  in  terms 
of  these  “essential  states”  the  Schrodinger  equation  be¬ 
comes 


b(t)=  ~2  Cna„(t)e‘  , 

n 

(2.1a) 

bj(t)=  — i(<o—a)0)bj(t )  — 2  Cnan(t)e'k''  %J  , 

n 

(2.1b) 

an(t)=  —i  (co„  —coQ)an(t)  +  Cnb(t)e  *k"  *° 

+  CnJtbj(t)e~iK'Xj  . 

(2.1c) 

j 


We  have  set  the  energy  scale  such  that  the  state  with  the 
atom  at  x0  excited  has  energy  zero.  The  atom-field  in¬ 
teraction  Hm{  has  been  taken  to  be  the  “electric-dipole 
form” 


(2.2) 


The  electric  field  has  been  expanded  in  a  complete  set  of 
free-space,  plane-wave  modes  with  associated  photon  an¬ 
nihilation  operators  an.  Actually,  for  our  purposes  it  is 
very  convenient  to  restrict  our  considerations  to  modes 
with  wave  vectors  k„  parallel  to  the  z  axis  of  Fig.  1.  Thus 
we  take 


Cn  =  ^r(2irHcon/AL)'/2  ,  (2.3) 

n 

where  A  is  an  effective  area,  L  is  the  length  along  the  z 
axis  of  our  quantization  box,  and  p  is  the  magnitude  of 
the  transition  dipole  moment  of  each  two-state  atom.  The 
operators  5 ,  and  dj  in  (2.2)  are,  respectively,  the  raising 


-J/ 
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Excited  atom  at  >T0 


FIG.  1.  Excited  atom  at  Xo  located  a  distance  d  from  a 
dielectric  half-space  of  N  atoms  per  unit  volume. 


Dielectric  medium 
of  N  atoms/volume 


and  lowering  operators  for  the  yth  two-state  atom.  For 
simplicity  we  take  each  field  mode  to  be  linearly  polarized 
along  the  direction  of  the  identical  dipole  moments. 

Equation  (2.1c)  can  be  used  to  formally  eliminate  the 
amplitudes  an(t)  from  (2.1a)  and  (2.1b).  The  result  is  the 
coupled  set  of  equations 


M/)=-2Cn  Sldt'bi 


t')e 


i  ia>„  -  cj.,11  /  -  f  i 


-2  2Cy^Z""Z',fl'd!'bJU')e“^ t 
i  J 


(2.4a) 


bj(t)  = 


—  Hco—a)  o  )bj  ( t ) 


OJ  —  11  t  -  I  I 


I  01  -01,,  II  I' -II 


(2.4b) 


Here  we  have  used  our  assumption  that  the  atom  at  x0  is 
excited  at  t  =0,  so  that  M0)=  1,  bj(0)=a„(0)  =  0. 

In  the  limit  in  which  the  length  L  of  our  quantization 
box  goes  to  infinity  we  have 


i(o>n  — <u0}(/'  — f)  (A 
?  = 


2v -ft 


AL 


i(«y_  —  o)nur  —  f) 


p  dkwe,>cj'-a°)u'-" 

UAL  2rr  J 


(2.5) 


where  |  k'  |  =<y'/c.  If  it  is  assumed  that  only  modes  with 
frequencies  interact  strongly  with  the  atoms,  one 

might  replace  the  integration  variable  k'  by  w'/c,  allow¬ 
ing  co'  to  take  on  negative  values.  This  is  a  good  approxi¬ 
mation  for  our  purposes.  Thus 


2c»V  v 


-ianU'-t) 


fiAc 


f  "  dh)'(o'eia,i,'-,) 

»  —  oc 


=  _JuLe-i"o "'-''A  p  da-eM-» 

fiAc  6 1'  •'-« 

fiAc  bt 


(2.6) 


In  this  approximation  we  have 

2 Cl 


2mn2 

fiAc 


f  dt'b{t')e~ia 0< 

J  o 


—  f) 


5  7TUO)0 

HaT™- 


(2.7) 


■ry.'yy. 
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We  are  ignoring  an  infinite  term  corresponding  to  a  theory,  specialized  to  the  case  of  field  modes  propagating 
single-atom  frequency  shift.  We  have  also  made  use  of  along  the  z  axis, 
the  fact  that  b(t)  is  slowly  varying  compared  with  e  0 .  Similarly 
Our  approach  is  just  a  variant  of  Weisskopf-Wigner 


—  —  f"  dk'o>cos(k'z)eUa,'-a°)[,'-t) 

AL  j  2xr  J  -<* 

/°C  I  ((O  —  0>qK/'  —  /) 

do)  u>  cos (to  z/c)e 

—  00 


Ha)'  —  —  t) 


— [6(/'-f  +z/c)  +  S(r'  — /  - z/c )] 
bt 


and  so 

zcyk-^'  f'dt'bu^-^''-0 


bj{t  —lj/c)Q[t  —lj/c) 


-  r^^eik"l'biit  —  I, /c )©( t -L/c)  ,  (2.9) 

fiAc  '  1  1 

where  0  is  the  unit  step  function,  k0=a>0/c,  and 

l j  —  i  Z()  Zj  i  . 

Combining  these  results  in  Eqs.  (2.4),  we  have  the 
delay-differential  equations" 

bU)=  -Kb(t)-K  2  e,k°^bj(t  -lj/c)Q(t  —  lj/c)  , 
i 

(2.10a) 

bj(t)=  —  i  (to  —  (t)0)bj(t) 

—  Ke'k",/b(t  ~lj/c)Q(t -Ij/cl-Kbjit) 
-K^e ,k"1'’ b,(t  - Ijj /c )©( r  - /c )  , 


(2.10b) 


s  — e ika>ib (t -21 . /cm t- 21. /c ) 

(O—COq 


and,  from  (2.10a) 


b(t)=-Kb(t)--^—yde1‘ku‘'bU-2Ij/c) 

co-to0 


XGU-2/y/c).  (2.13) 

This  delay  differential  equation  describes  the  effect  of  the 
dielectric  in  Fig.  1  on  the  probability  amplitude  b  for  the 
atom  at  x0  to  be  excited. 

Now  we  pass  to  the  limit  in  which  the  dielectric  con¬ 
tains  NAdz  atoms  in  the  slice  [z,z +dz],  making  the  re¬ 
placement 

2  e  2ik°ljb  ( t  —  21  j  /c  )©(  t  —  21  j  /c ) 
i 

k,  a  r  * 


where  /,,  =  z, and  K  =  TTf2:to()/^Ac.  Equations  -+NA  f  *  dz'e2‘knl*  ~z°' 

(2.10)  express  the  effects  of  atom-atom  coupling  in  a  way 

that  displays  explicitly  the  retarded  nature  of  the  elec-  xM(  -  2(z'-zo)/c)0(i  -  2(z'-z0)/c)  . 

tromagnetic  interaction. 

Since  we  are  interested  in  the  case  in  which  only  the  (2.14) 

atom  at  x, ,  is  initially  excited,  it  is  reasonable  to  suppose 

that  the  coupling  b,~bj  of  probability  amplitudes  for  Partial  integration  yields  for  the  integral  the  approximate 
atoms  inside  the  dielectric  is  small  compared  with  the  expression 

b  bt  coupling  involving  the  atom  at  x0.  If 

•o  -o,,  »K.  furthermore,  we  may  write  _  — !—  el,k"db(t -2d/c)Q(l -2d /c)  ,  (2.15) 


b^n-  -Ke‘k"lj  f  dt'hU'  -Ij/cmt' -lj/c) 


'  (O  -  <iA,"  t'  1 1 


^  - Ij/cmt  —lj/c)  .  (2.11)  with  (2.13),  we  obtain 

10  -  t:),t 


--l—e2'k"dbU  —  2d /c)Q{ t  —2d/c)  ,  (2.15) 

2  ike 

where  d  =1  — z0  is  the  distance  of  the  atom  at  x()  from  the 
dielectric  interface.  The  approximation  (2.15)  uses  again 
the  fact  that  bit)  varies  slowly  compared  with  e  ,0>0 ,  and 
we  ignore  rapidly  oscillating  terms.  Combining  this  result 


This  "adiabatic  following"  approximation  is  obtained  by 
partial  integration,  making  use  of  the  fact  that  b(t)  is 
slowly  varying  compared  with  e  "  ,  because  of  our 
assumption  that  K  <r<r  <o  —  tolt  . 


bU)=-Kb(t)+K  HhOLfi*  em"db(t-2d/c) 

(0  —  fU() 


xO(t~2d/c)  . 


□□ 
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The  refractive  index  of  a  dielectric  of  N  two-state 
atoms  per  unit  volume,  each  of  transition  frequency  a>  and 
transition  dipole  moment  /(,  is  given  by 


n  («()!'  - 


4  TT.\fl  :io/fl 

ft)"  —  a>o 


(2.17) 


if  local  field  effects  are  ignored,  as  is  permissible  for  a  di¬ 
lute  medium.  Of  course,  this  is  consistent  with  our  earlier 
neglect  of  the  last  term  in  (2.10b),  which  corresponds  to 
atom-atom  interactions  within  the  dielectric.  For  io  =  a>0 
and  it  uj,,)  =  1,  therefore,  we  have 


.Vrjr/2#  ^  n  -  1 

<:J  —  W,  |  II  +  1 


(2.18) 


and 

bills  —K[biO->-Ri‘~'k,'Jbit  —  2d  /c)Q(t  —2d/c)]  , 

(2.19) 

where  R  =  —  in  —  D/in  + 1 )  is  the  reflection  coefficient 
according  to  the  Fresnel  formula  for  normal  incidence. 

The  solution  of  (2.19)  with  the  initial  condition  b(0)=  1 
is 

(  -KRe2‘k"din ,  „  ,  ,  „ 

bit)  =  T - : - it  -2 nd/e) 

.  n ! 


Xe-Kit-lnd/c,QU  _2 nd/c)  .  (2.20) 

Similar  solutions  were  discussed  some  time  ago  for  the 
resonant  two-atom  interaction.12  When  2d /c -+0  (2.20) 
reduces  to 

b(n  =  e-K‘'  +  Re'‘k"  '• ,  (2.21) 

which  also  could  have  been  deduced  from  (2.19).  From 
this  result  we  can  see  that  K[  1  +/?  cos(2k0<f)J  and 
KR  sinilkod)  represent  a  decay  rate  and  frequency  shift, 
respectively,  of  the  atom  near  the  dielectric,  although 
these  expressions  are  unrealistic  in  the  sense  that  we  have 
only  included  modes  propagating  parallel  to  the  z  axis  in 
our  analysis.  The  more  realistic  expressions  including  all 
field  modes  are  easily  derived.1’  Our  point  here  is  that 
(2.20)  describes  how  the  atom-dielectric  interaction  builds 


up,  as  it  were,  by  “virtual  photon  exchanges.”  The  re¬ 
striction  to  modes  propagating  parallel  to  the  z  axis  be¬ 
comes  somewhat  more  meaningful  in  Sec.  III. 

Note  that  in  effect  a  layer  of  atoms  of  depth 
sc  k u  1  sc A.0  gives  rise  to  the  reflection  coefficient.  That  is, 
the  reflection  coefficient  arises  mainly  from  atoms  near 
the  surface  of  the  dielectric  interface.  This  is  seen  from 
the  approximation  (2.15)  to  the  integral  appearing  in 
(2.14).  From  a  rigorous  point  of  view,  however,  the  re¬ 
flection  coefficient  has  contributions  from  all  the  atoms 
comprising  the  dielectric.  We  discuss  this  point  further  in 
Sec.  IV. 

III.  ATOM  IN  A  FABRY-PEROT  RESONATOR 

We  first  consider  an  atom  between  two  perfectly  con¬ 
ducting  walls  (Fig.  2).  Then  we  will  generalize  to  the  case 
of  mirror  reflectivities  R^t-  —  1.  We  could  proceed  as  in 
the  preceding  section,  but  instead  we  will  follow  a  slightly 
different  approach,  expanding  the  field  in  mode  functions 
sin(k„z),  k„=nn/L.  The  Schrodinger  equation  now 
takes  the  form 

bit)  —  —  D  ]£ft)i/26„(/)sin(A'„z0)  ,  (3.1a) 

n 

b„it)  =  —  ii(o„  —(o0)bnit)  +  Da>]n/2bit)smiknzQ)  ,  (3.1b) 

where  D  =(4n^2 /fiAL)xn  and  again  we  have  employed 
the  essential-states  approximation,  bit)  being  the  ampli¬ 
tude  for  the  state  in  which  the  atom  is  excited  and  there 
are  no  photons  in  the  field,  and  bnit)  the  amplitude  for 
one  photon  in  mode  n  and  the  atom  deexcited.  Using  the 
formal  solution  of  the  second  equation  in  the  first,  we 
have  [/)„(0)=0,  M0)  =  1] 

bit)—  —D2  J  dt’bit’)  2  w,  sin2(k„z0) 

n  =0 


where  k„—con/c,  and  A  =  7 tc/L  is  the  mode  (an¬ 

gular)  frequency  separation.  Now  by  the  same  sort  of  ar¬ 
gument  used  in  Sec.  II  we  let  the  ton  take  on  negative 
values,  arguing  that  frequencies  far  removed  from  the 
atomic  resonance  frequency  cannot  make  a  significant 
contribution.  Thus  we  replace  (3.2)  by 


hitt  iD:  f  dt'b 

J  I) 


*  2  sin 2lknz,)e‘^r  “ 

n  --  x 


I  r,  ’  f  '  I  ,  ,  "''n  I'  f> 

£>-  dt  bit  ie 

y  j  o 


bt' 


t'  1 1  J_  In ^  1  1  4  ^zn^c>  1  {  — 


(3.3) 


From  the  Poisson  summation  formula 


2  e‘nx  =  j  tox/2 n-n)  . 


(3.4) 


therefore,  we  obtain  after  some  straightforward  manipulations  as  above  the  delay-differential  equation 
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bit )=  -Ul:T  jftUH-  2  e,nU{)’b(t -nT)Q(t -nT)- \ehk^  2  ena°rb(t  —  T0  —  nT)QU  —  T0  —  nT) 

n  =  1  n  =  0 

-±e  llk"!"  2  e‘no,"'b(t  +  T0-nT)Oit  +  T0-nT)  ,  (3.5) 

n  —  1 

where  Vl:  =  D2<j0  =  4-^i'n>„/fiAL,  k0  =  ioit/c,  T  =  2L/c,  and  Tu  =  2z0/c.  This  equation  displays  retardation  effects  due 

to  the  presence  of  two  walls,  and  as  such  has  a  much  more  complicated  delay-time  structure  than  (2.19). 

Writing  +  where  «  is  the  frequency  of  the  cavity  mode  closest  to  the  atomic  frequency  to0,  we  have 

iiotiV  i  to  +  A,,1  ZL  /c  i A,,/  ,  , 

e  =e  '  =e  (3.6) 

as  a  result  of  the  mode  condition  sinWjL  /c)  =  0.  Therefore  we  may  replace  e‘nu>{)1  by  e'"^'1  in  (3.5). 

To  account  for  mirror  reflectivities  R  <  1  we  modify  (3.5)  based  on  the  results  of  Sec.  II.  Assuming  both  mirrors 
have  the  same  reflectivity  R  =  —  in  —  1 )/(«  +  1 ),  and  associating  with  each  “bounce"  off  a  mirror  a  factor  R,  therefore, 
we  replace  (3.5)  by 

i>it)  =  -yll-Y  t  6(/>+  2  (R:c‘S"'  )nbit -nT)6(t -nT)  +  jRehk"Z"  2  )nb  It  -  T0 -nT)0U  -  T0-nT) 


+  7 R-'e  2'Vo  2  iR2e‘^')nb(t +T0-nT)6(t +T0-nT)  .  (3.7) 

n  =  1 

Equation  (3.7)  is  quite  general  in  that  it  includes  possible  effects  of  all  longitudinal  modes,  as  well  as  mirror  reflectivities 
/?7=  —  1.  In  general,  however,  the  time  dependence  prescribed  by  (3.7)  is  rather  complicated.  For  this  reason  it  is 
worthwhile  to  focus  our  attention  on  some  special  cases,  and  show  that  some  well-known  results  can  be  recovered  from 

(3.7). 

A.  Rabi  oscillations 

For  perfectly  reflecting  mirrors,  R  —  —  1.  For  (IT  «  1  and  t»T  we  have  in  this  case 
6U)=-4nVA>'  T  2  e~ilko"~nT)b(t  -nT)Q(t  -nT)-\elik°t0T  2  e~iKP~nT'b(t -T0-nT)Q(t -T0-nT) 

n  =  I  (1=0 


A 


—  ye  2ik°zor  2  e  'A,,t  nT'b(t  +  T0  —  nT)O(t  +T0  —  nT) 

n  =  1 

-  -  4fiVV(  1  -  je2'Vfl-  4e  ~2'Vo)  /'  dt'b(f)e'^' 

=  — ^2sln2(koZole,An,  fQ  dt'b(t')e  ,A 11 ' 


bit)  -  i  &,£(!)  +  —  6(f)  =  0 
4 

with  /,  =  211  sin(  kuzt) ).  Thus 

tb ( / >  =  e ' A"f  I  cos[  4- ( Ao -+-  A2)l/2r] 


TTT  73  sin[  T ( Aft  +  >  <] 

r  t  -  \  l '  - 


(A«  +  a 


photon  bounce  frequency  T~]  =c /2L  is  much  greater 
than  the  Rabi  frequency.  Note  that  when  this  condition  is 
satisfied  we  also  have  n27'«7'~1  =  A,  i.e.,  the  cavity¬ 
mode  spacing  is  large  compared  with  the  spontaneous  de¬ 
cay  rate  and  therefore  the  natural  linewidth.14  Further¬ 
more  n-r «n  means  that  the  spontaneous  emission  rate 
is  small  compared  with  the  Rabi  frequency.  Thus  it  is  not 
surprising  that  Rabi  oscillations  occur  in  this  limit.  This 


which  displays  the  well-known  Rabi  oscillations  for  an 
atom  interacting  with  a  single  field  mode.  (In  this  case 
they  may  be  termed  “vacuum"  Rabi  oscillations.)  For  ex¬ 
act  resonance,  Ao  =  0,  we  have  bit)  =  cos(  jkt),  and  the 
population  difference 

hit) ;  2  —  (1  —  bit) , 2 )  =  cos( Ar )  .  (3.11) 

Thus  our  general  delay-differential  equation  (3.7) 
predicts  Rabi  oscillations  when  (IT  «\,  i.e.,  when  the 


FIG.  2.  Excited  atom  between  two  dielectric  walls,  located  a 
distance  z„  front  the  nearest  wall. 
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single-mode  limit  is  the  same  limit  for  which  retardation 
effects  are  negligible. 

B.  Damped  single-mode  cavity 

Now  let  us  consider  the  same  single-mode  limit  for 
R  x  —  1.  Writing 

:  .A „/  .iA „*r  1  ln«  -  ■  I 

K  e  =e 


with  y  =  —  ( c  /2L  >ln/? :  >0  the  field  damping  rate  due  to 
imperfectly  reflecting  walls,  we  can  perform  the  same  ma¬ 
nipulations  that  led  to  (3.9)  to  obtain 

bit)  —  /(A(>+  iy)bit) -f  —b (f)s0  .  (3.13) 

4 

For  Ao  =  0  and  k  »}',  for  instance,  we  have 

bit)  :-re  _y'cos~(  -rAf)  .  (3.14) 

For  y  >■>/,,  on  the  other  hand, 

bit)  -  .  (3.15) 

These  results  agree  with  those  of  Purcell*  and  Sachdev.fl 
In  fact,  the  solution  of  (3.13)  with  A0  =  0  is  exactly 
equivalent  to  Sachdev’s  general  solution  (5.2)  obtained  by 
a  different  approach.  We  therefore  refer  the  reader  to 
Ref.  9  for  a  discussion  of  the  single-mode  case  with  cavity 
damping.  (See  also  Ref.  10.) 

C.  Damped  multimode  cavity 

If  T  is  increased  the  solution  of  (3.7)  has  a  complicated 
delay-time  structure,  as  discussed  earlier  for  a  similar  but 
somewhat  simpler  quantum  delay-differential  system.1* 
In  particular,  if  T  is  large  then  the  Rabi  frequency  is  large 
compared  with  the  photon  bounce  frequency.  The  atom 
can  emit  spontaneously  before  feeling  the  presence  of  the 
cavity,  and  later  it  can  absorb  the  emitted  photon,  reemit, 
etc.,  without  any  coherent  Rabi  oscillations  associated 
with  the  single-mode  limit  517' «1.  Since  the  solutions 
in  this  case  resemble  those  shown  in  Ref.  15,  we  will  not 
take  the  time  here  to  display  them  graphically. 


IV.  REMARKS  ON  THE  EWALD-OSEEN 
EXTINCTION  THEOREM 

According  to  the  classical  Ewald-Oseen  extinction 
theorem,16  the  polarization  induced  in  a  dielectric  medium 
produces  in  the  medium  a  Field  that  consists  of  two  parts. 
One  part  exactly  cancels  the  incident  Field  inside  the 
medium,  whereas  the  other  propagates  inside  the  medium 
at  the  phase  velocity  c/n.  The  Field  radiated  out  from  the 
medium  is  just  the  reflected  Field,  with  amplitude  deter¬ 
mined  by  the  Fresnel  reflection  coefFicient. 

In  Sec.  II  we  obtained,  starting  from  the  Schrodinger 
equation,  the  correct  Fresnel  reflection  coefFicient  for  noi- 
mal  incidence,  assuming  nsl.  The  approximation  (2.15) 
to  the  integral  appearing  in  (2.14)  indicates  that  a  number 
of  atoms  xzNA  /k0z=NAk0  contributes  to  the  reflected 
Field.  That  is,  atoms  in  a  layer  of  depth  s=A.0  at  the  sur¬ 
face  of  the  dielectric  cooperate  to  produce  the  reflected 
Field.  This  result  is  consistent  with  classical  argu¬ 
ments.1718 

To  describe  the  extinction  of  the  incident  Field  inside 
the  medium,  consider  the  steady-state  solution  oF  Equa¬ 
tion  (2.10b), 

bj  =  -(iK/^0)beiko>t-(iK/^o)1bieik°l‘J  .  (4.1) 

i*i 

Here  we  have  again  used  the  nonresonance  assumption 
j  A0|  =  |  w0-w|  »K  to  replace  A0-f  iK  by  A0.  The 
steady-state  assumption  is  useful  in  order  to  focus  our  at¬ 
tention  on  a  single-frequency  component  at  a  time  (as  in 
the  classical  Ewald-Oseen  extinction  theorem).  If  the 
second  term  on  the  right  side  of  (4.1)  were  absent,  we 
would  have 

i  bj  1 2  =  (/f2/Ao)  |  A  | 2  (4.2) 

for  the  probability  of  exciting  some  atom  j  in  the  medium. 
We  now  ask  how  this  result  is  modified  by  the  presence  of 
the  last  term  in  (4.1),  In  other  words,  what  is  the  proba¬ 
bility  that  the  initially  excited  atom  can  excite  an  atom 
within  an  entire  dielectric  medium  of  atoms? 

In  the  continuum  limit  we  replace  the  summation  in 
(4.1)  by  an  integral  over  z,  as  in  Sec.  II.  Writing  biz)  in¬ 
stead  of  bj ,  we  then  have  the  following  integral  equation 
for  biz): 


biz)—  -  UK /A)belkl  —  HK/A„).\A  f  dz'biz')e'k  z~z 

J  0 

=  -HK/\(,)beikz-iiK/At))SA  [  f'dz'b(z’)e'k'z-z  '+  dz'biz')eik'z’ ~z' 


(4.3) 


where  we  have  written  k  for  kt).  To  solve  this  equation 
we  seek  a  solution  of  the  form 

biz)  =  Celk  z  .  14.4) 

where  C  and  k'  are  constants  to  be  determined.  Using 
this  form  in  (4.3),  and  performing  the  integrals,  we  obtain 

Ce‘kz=  -UK/b„)he,kz 

e,kz-e,kz  _e 
k'  —  k  k'  4- k 


and,  equating  coefFicients  of  elkz  and  c‘k  z. 

k'2-k2=  -2kKNA  /A„  , 

(4.6a) 

C  =iik'-k)/NA  . 

(4.6b) 

Equation  (4.6a)  gives  k'  =  kn,  where 

n 2  -  1  =  —  IKS  A  /k  A0 

(4.7a) 

or 


( A.V.4C/A,,) 


(4.5) 


n  s  1  +nfi2N /fHw  —  coo) 


(4.7b) 
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for  n  —  1  «  1.  n  is  just  the  refractive  index  for  light  of 
frequency  <j0  in  a  medium  of  .V  two-state  atoms  per  unit 
volume,  each  with  transition  dipole  moment  p  and  transi¬ 
tion  frequency  to.  [Equation  (2.18).] 

Equation  (4.6b)  is  the  condition  that  the  incident  field 
in  the  medium  is  exactly  cancelled  by  the  part  of  the  di¬ 
pole  field  in  the  medium  that  varies  as  exp(/Az).  [See  Eq. 
(4.5'.]  Using  our  result  for  k',  (4.6b)  gives 

C  =  — 22—  (  -  iKb  / A(1 )  s  -  -2'L- (  - iKb  / A0 )  .  (4.8) 

n  +  1  n + 1 

From  Eqs.  (4.2)  and  (4.5)  we  see  that  the  amplitude  to  ex¬ 
cite  any  atom  in  the  medium  is  a  factor  2 n  /( n  +  1 )  times 
the  amplitude  to  excite  that  atom  if  it  were  in  free  space. 
This  factor  is  just  the  Fresnel  transmission  coefficient  (for 
normal  incidence,  because  in  our  model  we  only  allow 
plane-wave  modes  propagating  normally  to  the  dielectric 
interface). 

These  results  provide  a  fully  microscopic ,  quantum- 
mechanical  basis  for  the  Ewald-Oseen  extinction  theorem, 
starting  from  the  Schrodinger  equation.  Equation  (4.3)  is 
a  quantum  analogue,  in  terms  of  probability  amplitudes, 
of  the  classical  superposition  principle  for  the  field. 


The  extinction  theorem  is  a  nonlocal  boundary  condi¬ 
tion  that  the  field  must  satisfy. IS  Physically,  the  cancella¬ 
tion  of  the  incident  field  is  often  regarded  as  “caused  by 
the  dipoles  on  the  boundary  of  the  medium,”20  since  in  the 
classical  macroscopic  approach  the  term  that  cancels  the 
incident  field  can  be  cast  in  the  form  of  a  surface  integral 
over  dipole  sources.  Our  approach  shows  how  all  the  di¬ 
pole  contributions  add  up  in  such  a  way  that  the  cancella¬ 
tion  is  effectively  due  to  dipoles  within  a  depth  approxi¬ 
mately  equal  to  k  at  the  surface.  In  particular,  “The  (re¬ 
flected)  radiation  comes  from  everywhere  in  the  interior, 
but  it  turns  out  that  the  total  effect  is  equivalent  to  a  re¬ 
flection  from  the  surface.”111 
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